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Abstract 

Natural objects can be subject to various transformations yet still preserve 
properties that we refer to as invariants. Here, we use definitions of affine in- 
variant arclength for surfaces in in order to extend the set of existing non-rigid 
shape analysis tools. In fact, we show that by re-defining the surface metric as 
its equi-affine version, the surface with its modified metric tensor can be treated 
as a canonical Euclidean object on which most classical Euclidean processing and 
analysis tools can be applied. The new definition of a metric is used to extend 
the fast marching method technique for computing geodesic distances on surfaces, 
where now, the distances are defined with respect to an affine invariant arclength. 
Applications of the proposed framework demonstrate its invariance, efficiency, and 
accuracy in shape analysis. 



1 Introduction 

Modeling 3D shapes as Riemannian manifold is a ubiquitous approach in many shape 
analysis applications. In particular, in the recent decade, shape descriptors based on 
geodesic distances induced by a Riemannian metric have become popular Notable 
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examples of such methods are the canonical forms Q and the Gromov-Hausdorff ||9l 
[131 |2I and the Gromov-Wasserstein |12| frameworks, used in shape comparison and 
correspondence problems. Such methods consider shapes as metric spaces endowed 
with a geodesic distance metric, and pose the problem of shape similarity as finding 
the minimum-distortion correspondence between the metrics. The advantage of the 
geodesic distances is their invariance to inelastic deformations (bendings) that preserve 
the Riemannian metric, which makes them especially appealing for non-rigid shape 
analysis. A particular setting of finding shape self-similarity can be used for intrinsic 
symmetry detection in non-rigid shapes ifTTl [23] fTTl l22l . 

The flexibility in the definition of the Riemannian metric allows extending the in- 
variance of the aforementioned shape analysis algorithms by constructing a geodesic 
metric that is also invariant to global transformations of the embedding space. A par- 
ticularly general and important class of such transformations are the ajfine transfor- 
mations, which play an important role in many applications in the analysis of images 
lfT4l and 3D shapes |8|. Many frameworks have been suggested to cope with the ac- 
tion of the affine group in a global manner, trying to undo the affine transformation in 
large parts of a shape or a picture. While the theory of affine invariance is known for 
many years [4| and used for curves ifTSl and flows |fT9l . no numerical constructions 
applicable to manifolds have been proposed. 

In this paper, we construct an (equi-)affine-invariant Riemannian geometry for 3D 
shapes. By defining an affine-invariant Riemannian metric, we can in turn define affine- 
invariant geodesies, which result in a metric space with a stronger class of invariance. 
This new metric allows us to develop efficient computational tools that handle non-rigid 
deformations as well as equi-affine transformations. We demonstrate the usefulness of 
our construction in a range of shape analysis applications, such as shape processing, 
construction of shape descriptors, correspondence, and symmetry detection. 



2 Background 

We model a shape {X, g) as a compact complete two-dimensional Riemannian man- 
ifold X with a metric tensor g. The metric g can be identified with an inner product 
(•, ■)x ■ TxX X TxX — > M on the tangent plane T^X at point x. We further assume 
that X is embedded into M.^ by means of a regular map x : [/ C ^ M^, so that the 
metric tensor can be expressed in coordinates as 

dui duj ' 

where Ui are the coordinates of U. 

The metric tensor relates infinitesimal displacements in the parametrization domain 
U to displacement on the manifold 

dp^ — giidui^ + 2gi2duidu2 + g22du2^ . (2) 

This, in turn, provides a way to measure length structures on the manifold. Given a 
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curve C : [0, T] — > X, its length can be expressed as 

m = f {c[t),c[t))'^i^dt, (3) 

where C denotes the velocity vector. 

2.1 Geodesies 

Minimal geodesies are the minimizers of ^(C), giving rise to the geodesic distances 

dxix,x')^ min 1(C) (4) 

where T{x,x') is the set of all admissible paths between the points x and x' on the 
surface X, where due to completeness assumption, the minimizer always exists. 

Structures expressible solely in terms of the metric tensor g are called intrinsic. 
For example, the geodesic can be expressed in this way. The importance of intrinsic 
structures stems from the fact that they are invariant under isometric transformations 
(bendings) of the shape. In an isometrically bent shape, the geodesic distances are 
preserved - a property allowing to use such structures as invariant shape descriptors 
0. 

2.2 Fast marching 

The geodesic distance dx{xo, x) can be obtained as the viscosity solution to the eikonal 
equation \\Vd\\2 ~ 1 with boundary condition at the source point d{xo) — 0. In the 
discrete setting, a family of algorithms for finding the viscosity solution of the dis- 
cretized eikonal equation by simulated wavefront propagation is called /aif marching 
methods 121^ Ji)J. On a discrete shape represented as a triangular mesh with N vertices, 
the general structure of fast marching closely resembles that of the classical Dijkstra's 
algorithm for shortest path computation in graphs, with the main difference in the up- 
date step. Unlike the graph case where shortest paths are restricted to pass through the 
graph edges, the continuous approximation allows paths passing anywhere in the mesh 
triangles. For that reason, the value of d{xQ,x) has to be computed from the values of 
the distance map at two other vertices forming a triangle with x. Computation of the 
distance map from a single source point has the complexity of 0{N log N). 

3 Affine-invariant geometry 

An ajfine transformation x M> Ax + b of the three-dimensional Euclidean space can 
be parametrized using twelve parameters: nine for the linear transformation A, and 
additional three, b, for a translation, which we will omit in the following discussion. 
Volume-preserving transformations, known as special or equi-affine are restricted by 
det A = 1. Such transformations involve only eleven parameters. 
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dC 
dp 


dui 
dp 


dp 


d^C 
dp"^ 


d'^ui 
dp"^ 


+ x ^^^^ 
dp'^ 



The equi-affine metric can be defined through the parametrization of a curve on the 
surface H] [3] |4l |6] [15] [I9j . Let C be a curve on X parametrized by p. By the chain 
rule, 

(5) 

/ dui \ ^ dui du2 ( du2 ^ ^ 

xii — — +2x12-5— -5 '"^22 -1— 

\ dp J dp dp \ dp 

where, for brevity, we denote Xi = ^ and x.^ — ^ . As volumes are pre- 
served under the equi-affine group of transformations, we define the invariant arclength 
p through 

/(X)dct(xi,X2,Cpp) = 1, (6) 

where f{X) is a normalization factor for parameterization invariance. Plugging (|6]l 
into (|6]l yields 

dp'^ = /(X) det(xi, X2, Xiidw^ + 2xi2(iui(iu2 + X22<iu2) 

= f{X) (giidul + 2gi2duiU2 + g22dul) . (7) 

where gij = det (xi , X2 , x^^ ) . 

In order to evaluate f{X) such that the quadratic form dTjl will also be parame- 
terization invariant, we introduce an arbitrary parameteriation ui and U2, for which 
Xi = and Xij ~ q^.q^. ■ The relation between the two sets of parameterizations 
can be expressed using the chain rule 

Xl — Xu-^ — Xu-^Uiui + X„2U2aj (8) 

X2 — Xu^ — Xu-i^Ulf^^ + Xu2U2u2- 

It can shown [4 1 using the Jacobian 



that 

giidul + 2gi2duidu2 + g22du\ = [gudu^ + 2gi2dudv + .922^"^) | J|, (10) 

and 311322 - 3?2 = {911922 - gh)"^, where = det(xi, X2, Xy). From ([TOll we 
evaluate the affine invariant parameter normalization f{X) — |g|~^^'', and define an 
equi-affine pre-metric tensor ||4l lT9i 

= 9ij\g\~^^'^- (11) 

The pre-metric tensor ( fTTT i defines a true metric only for strictly convex surfaces H; 
A simiak problem appeared in equi-affine curve evolution where the flow direction was 
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Figure 1 : The three neighboring triangles together with the central one are unfolded flat 
to the plane. The central triangle is canonized into a right isosceles triangle while the 
rest of its three neighboring triangles follow the same planar affine transformation. Fi- 
nally, the six surface coordinate values at the vertices are used to interpolate a quadratic 
surface patch from which the metric tensor is computed. 




Figure 2: Geodesic level sets of the distance function computed from the tip of the tail, 
using the standard (left) and the proposed equi-affine (right) geodesic metrics. 

determined by the curvature vector In two dimentions we can encounter non-positive 
definite metrics in concave, and hyperbolic points. We propose fixing the metric by 
flipping the main axes of the operator, if needed. In practice we restrict the eigenvalues 
of the tensor to be positive, and re-avaluate it. From the eigendecomposition G = 
UrU^ of g in matrix notation, where U is orthonormal and F = diag{7i,72}, we 
compose a new metric G, such that G = U|rjU^ is positive definite and equi-affine 
invariant. 



5 



Standard Equi-affine 




Standard Equi-affine 




Figure 3: Distance maps from different source points calculated using the standard 
(second to fourth columns) and the proposed equi-affine geodesic metric (fifth to sev- 
enth columns) on a reference surface (first and third rows) and its affine (second row) 
and isometric deformation-naffine transformation (fourth row). Thirds and sixth rows 
show the global histogram of geodesic distances before and after the transformation 
(green and blue curves). The overlap between the histograms is an evidence of invari- 
ance. 

4 Discretization 

We model the surface X as a triangular mesh, and construct three coordinate functions 
x(m, u), y{u, v), and z{u, v) for each triangle. While this can be done practically in 
any representation, we use the fact that a triangle and its three adjacent neighbors, can 
be unfolded to the plane, and produce a parameter domain. The coordinates of this 
planar representation are used as the parametrization with respect to which the first 
fundamental form coefficients are computed at the barycenter of the simplex (Figure 
[I)- Using the six base functions 1, u, v, uv, and we can construct a second order 
polynom for each coordinate function. This step is performed for every triangle of the 
mesh. 
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Calculating geodesic distances was well studied in past decades. Several fast and 
accurate numerical schemes ifTOl l20l l24l can be used off-the-shelf for this purpose. 
We use FMM technique, after locally rescaling each edge according to the equi-affine 
metric. 

The (affine invariant) length of each edge is defined by L^{dx, dy) — gudx^ + 
2gi2dxdy+g22dy^ ■ Specifically, for our canonical triangle with vertices at (0, 0), (1, 0) 
and (0, 1) we have 
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Each edge may appear in more than one triangle. We found that the average length is 
a good approximation, assuming the triangle inequality holds. In figures |2] and [3] we 
compare between geodesic distances induced by the standard and our affine-invariant 
metric. 



5 Results 




Figure 4: Voronoi cells generated by a fixed set of 20 points on a shape undergoing 
an equi-affine transformation. The standard geodesic metric (left) and its equi-affine 
counterpart (right) were used. Note that in the latter case the tessellation commutes 
with the transformation. 

The equi-affine metric can be used in many existing methods that process geodesic 
distances. In what follows, we show several examples for embedding the new metric in 
known applications such as voronoi tessellation, canonical forms, non-rigid matching 
and symmetry detection. 
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5.1 Voronoi tessellation 



Voronoi tessellation is a partitioning of {X, g) into disjoint open sets called Voronoi 
cells. A set of k points [xi G X),^^^ on the surface define the Voronoi cells {Vi)^^^ 
such that the i-th cell contains all points on X closer to Xi than to any other Xj in 
the sense of the metric g. Voronoi tessellations created with the equi-affine metric 
commute with equi-affine transformations as visualized in Figure|4] 

5.2 Canonical forms 

Methods considering shapes as metric spaces with some intrinsic (e.g. geodesic) dis- 
tance metric is an important class of approaches in shape analysis. Geodesic distances 
are particularly appealing due to their invariance to inelastic deformations that preserve 
the Riemannian metric. 

Elad and Kimmel Q proposed a shape recognition algorithm based on embedding 
the metric structure of a shape (X, dx) into a low-dimensional Euclidean spaces. Such 
a representation, referred to as canonical form, reduces the number of degress of free- 
dom by tanslating all deformations into a much simple Euclidean isometry group. 

Given a shape sampled at N points and an N x N matrix of pairwise geodesic 
distances, the computation of the canonical form consists of finding a configuration of 
N points zi , . . . , zat in such that 1 1 — Zj 1 1 2 ~ dx{xi, Xj ) . This problem is known 
as multidimensional scaling (MDS) and can be posed as a non-convex least-squares 
optimization problem of the form 

{zi,...,zn} = a,rgmm^\\\zi - Zj\\2 ~ dx{xi,Xj)\'^. (13) 

The invariance of the canonical form to shape transformations depends on the 
choice of the distance metric dx- Figure |5] shows an example of a canonical form 
of the human shape undergoing different bendings and affine transformations of vary- 
ing strength. The canonical form was computed using the geodesic and the proposed 
equi-affine distance metric. One can clearly see the nearly perfect invariance of the 
latter Such a strong invariance allows to compute correspondence of full shapes under 
a combination of inelastic bendings and affine transformations. 

5.3 Non rigid matching 

Two non-rigid shapes X, Y can be considered similar if there exists an isometric cor- 
respondence C C X X Y between them, such that Vx G X there exists y & Y with 
(x, y) ^ C and vice-versa, and dx{x, x') = dyiy, y') for all [x, y), (x' , y') G C, where 
dx , dy are geodesic distance metrics on X, Y. In practice, no shapes are truly iso- 
metric, and such a correspondence rarely exists; however, one can attempt finding a 
correspondence minimizing the metric distortion, 

dis(C) = max \dx{x, x') - dyiy^y')]- (14) 

(x,y)eC 
{x',y')eC 



8 



t 



i 



/ 
i 
h 

I 



t f / 



t 
\ 



r 
r 



X 

> 



T t :t 

X X X 

T y X 

X Xr ^ 



Figure 5: Embedding into of a human shape and its equi-affine transformations 
of varying strength. Classical scaling was used with a matrix of geodesic (left) and 
equi-affine geodesic (right) distances. In the latter case, canonical forms remain ap- 
proximately invariant up to a rigid transformation. 

The smallest achievable value of the distortion is called the Gromov-Hausdorff distance 
in between the metric spaces (X, dx) and (y, dy). 



and can be used as a criterion of shape similarity. 

The choice of the distance metrics dx , dy defines the invariance class of this sim- 
ilarity criterion. Using geodesic distances, the similarity is invariant to inelastic de- 
formations. Here, we use geodesic distances induced by our equi-affine Riemannian 
metric tensor, which gives additional invariance to affine transformations of the shape. 

Bronstein et al. |2| showed how ( fTSl l can be efficiently approximated using an 
optimization algorithm in the spirit of multidimensional scaling (MDS), referred to as 
generalized MDS (GMDS). Since the input of this numeric framework are geodesic 
distances between mesh points, all is needed to obtain an equi-affine GMDS is one 
additional step where we substitute the geodesic distances with their equi-affine equiv- 
alents. Figure |6] shows the correspondences obtained between an equi-affine trans- 
formation of a shape using the standard and the equi-affine-invariant versions of the 
geodesic metric. 

5.4 Intrinsic symmetry 

Raviv et al. IfTTl introduced the notion of intrinsic symmetries for non-rigid shapes 
as self-isometries of a shape with respect to a deformation-invariant (e.g. geodesic) 
distance metric. These self-isometries can be detected by trying to identify local min- 
imizers of the metric distortion or other methods proposed in follow-up publications 



Here, we adopt the framework of iflTl for equi-affine intrinsic symmetry detection. 
Such symmetries play an important role in paleontological applications IS] . Equi-affine 



dGYiiX.Y) = -infdis(C), 



(15) 



|[T6l|23l[ITl|22l. 
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Figure 6: The GMDS framework is used to calculate correspondences between a shape 
and its isometry (left) and isometry followed by an equi-affine transformation (right). 
Matches between shapes are depicted as identically colored Voronoi cells. Standard 
distance (first row) and its equi-affine-invariant counterpart (second row) are used as 
the metric structure in the GMDS algorithm. Inaccuracies obtained in the first case are 
especially visible in the legs and arms. 
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Figure 7; As the affine transformation becomes stronger, the quality of the symmetry 
detection decreases (B) when the standard geodesic metric is used. On the other hand, 
detection quality is hardly affected (C) by the transformations when using the equi- 
affine geodesic metric. 

intrinsic symmetries are detected as local minima of the distortion, where the equi- 
affine geodesic distance metric is used. Figure|2]shows that using the traditional metric 
we face a decrease in accuracy of symmetry detection as the affine transformation 
becomes stronger (the accuracy is the average geodesic mismatch with relation to the 
ground-truth symmetry). Such a decrease does not occur using the equi-affine metric. 

6 Conclusions 

We introduced a numerical machinery for computing equi-affine-invariant geodesic 
distances. The proposed tools were applied to applications, like symmetry detection, 
finding correspondence, and canonization for efficient shape matching. We have ex- 
tended the ability to analyze approximately isometric objects, like articulated objects, 
by treating affine deformations as well as non-rigid ones. As our analysis is based on 
local geometric structures, the affine group could in fact act locally and vary smoothly 
in space as long as it is encapsulated within our approximation framework. We intend 
to explore this direction and further enrich the set of transformations one can handle 
within the scope of metric geometry in the future. 
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